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Abstract 

We generalize the Callan-Harvey mechanism to the case of actions with 
a non local mass term for the fermions. Using a 2 + 1-dimensional 
model as a concrete example, we show that both the existence and 
properties of localized zero modes can also be consistently studied 
when the mass is non local. We derive some general properties from a 
study of the resulting integral equations, and consider their realization 
in a concrete example. 
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The Callan-Harvey mechanism PQ explains the existence and properties 
of fermionic zero modes that appear whenever the mass of a Dirac field 
in 2k + 1 dimensions (k = 1,2,- ■ ■) has a domain- wall like defect. The zero 
modes due to this phenomenon are localized, concentrated around the domain 
wall, and chiral from the point of view of the domain-wall world-volume (a 
2A;-dimensional theory). 

This mechanism has found many interesting applications, both to phe- 
nomenological issues [2] and theoretical elaborations E] • Dynamical do- 
main walls have been considered in [Sj, and the case of a (dynamical) super- 
symmetric model with this kind of configuration has been discussed in jB]. 
Remarkably, lattice versions of the zero modes, the so-called 'domain wall 
fermions' [7] have also been precursors of the overlap Dirac operator [S], a 
sensible way to put chiral fermions on a spacetime lattice. 

In spite of the fact that the Callan-Harvey mechanism has been extended 
in many directions, the existence of chiral zero modes has been so far only 
studied for the case of local mass terms, namely, those where the mass is 
a local function of the spacetime coordinates. In this article, we relax that 
condition and consider the case of a non local mass in 2 + 1 dimensions. This 
is a non trivial modification which does, however, arise naturally in some 
applications. For example, when the fermionic action is dressed by quantum 
effects, the corrections usually adopt the form of non local form-factors. 

Let us first review the standard Callan-Harvey mechanism Q when one 
considers a fermion field in 2+1 dimensions coupled to a domain wall like 
defect, with the (Euclidean) action 



We use Euclidean coordinates x = (xo, X\, x?) (xo is the Euclidean time) and 
$ = Ipd/i ; where the 7-matrices are chosen according to the convention: 



The local mass m(x) contains a topological defect; in the simplest case of a 
rectilinear static defect they have the characteristic shape: 



where o~(x-2) = sign(x 2 ). Therefore the domain wall, which is the interface 
between two regions with different signs for m(x), is the x± axis. 




(1) 



7o = 03 71 = a x 72 = cr 2 • 



(2) 




(3) 
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From the Dirac operator 

V = ft x + m(x 2 ) , (4) 

we can construct the hermitian operator 7i = T>^ T>. The form of 7i suggests 
the introduction of the adjoint operators 

a = d 2 + m(x 2 ) , a) = —d 2 + m(x 2 ) (5) 

in terms of which 

n = {a) a - $ 2 )P L + (aa f - $ 2 )P R (6) 

and 

V = (a+ 0)P L + (a f + 0)P R (7) 

where Pl = |(1 + 72), Pr = |(1 — 72). Expanding tf)(x) in the complete set 
of eigenf unctions of a^a and aa\ there appears a massless left fermion, 
localized over the domain wall. Its x 2 dependence is dictated by the fact 
that it is a zero mode of the a operator, and it dominates the low energy 
dynamics of the system. 

We want to generalize this phenomenon to include a non local mass. 
Namely, rather than (@J), the Dirac operator shall be 

V{x,y) = @ x 5(x-y) + M(x,y) . (8) 

Little can be said about the existence of a fermionic zero mode before we 
make some hypotheses to restrict the form of M(x,y). We assume that the 
system has translation invariance in the coordinates ) and that 

M(x, y) consists of a local domain wall like part, plus a non local term with 
a strength controlled by a parameter A: 

M(x,y) = m(x 2 ) 5(x - y) - X J e ik <*-®>y k (x 2 ,y 2 ) . (9) 

We are looking for a zero mode so that: 

(x\V\V) = J d 3 y t>(x,y)V(y) = . (10) 

Taking advantage of the translation invariance in x, we use 'separation of 
variables' to look for solutions of the form: 

*(x) = x(z)^(x 2 ) , (11) 
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where xi x ) is a massless spinor, which is left-handed from the point of view 
of the two-dimensional world defined by x, i.e., 

fa(x) = P rX (x) = . (12) 

There is an essential difference regarding the space of solutions to the equa- 
tions above in Euclidean and Minkowski spacetimes. Indeed, in Euclidean 
spacetime, it leads to analytic functions of Xq + ixx, while in the Minkowski 
case one has 'left-mover' solutions. Keeping this distinction in mind, we 
continue working with the Euclidean version. 

Substituting (|11|) into (|10|) and comparing with (J7J), we arrive to a non 
local version of the kernel for the annihilation operator 

a{x 2 ,y 2 ) = [d 2 + m(x 2 )] 5{x 2 - y 2 ) - A 7^(^2,2/2) • (13) 
The zero mode i{j(x 2 ), must then satisfy the equation 



(x\a\i(>) = [d 2 + m{x 2 )\ ip(x 2 ) - A / dy 2 >y k (x 2 , y 2 )ip(y 2 ) = . (14) 

J —00 

Following the method of variation of parameters, we use the ansatz 

ip(x 2 ) = ipo(x 2 )(p(x 2 ) (15) 
where ipo is the zero mode for the local part in ()14j). which satisfies 

[d 2 + m(x 2 )] ip (x 2 ) = i>o(x 2 ) = N exp[- / dsm(s)] , (16) 



and is a normalization constant. 

The function (p(x 2 ), which modulates ^'0(^2), must then satisfy the equa- 
tion 

/+00 
dy2 [ipo{x2)~ 1 1k(x 2 , y 2 ) ^0(2/2)] <p(y2) = . (17) 
-00 

By integrating over X2, the previous equation can be written in a more con- 
venient form as an integral equation 



<p(x 2 ) -A / dy 2 -f k (x 2 ,y 2 )ip(y 2 ) = (f(0) (18) 
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where the kernel 7^ is 

lk(x 2 ,y 2 ) = / dz 2 ip (z 2 y 1 >y k (z 2 ,y2) ip (y 2 ) ■ (19) 
Jo 

This is a homogeneous integral equation, which can be conveniently rewritten 
as an equivalent non-homogeneous set of equations. Indeed, introducing 
linear operators, (j!8|) can be rewritten as follows: 

(I-XT)tp = c (20) 

where 

(I<p)(x 2 ) = <p(x 2 ) , c=^(0), (21) 



(T(p)(x 2 ) = J dy 2 'j k (x 2 ,y 2 )(p(y 2 ) . (22) 

We see that the problem can be discussed in terms of (|2L)|). which is an 
inhomogeneous system of the Fredholm type 0. The condition c = (f(0) has 
to be verified, of course, after solving (|2*0|) for arbitrary c. 

The reason for this procedure is that, had we used the original homoge- 
neous system, we should have had to introduce non compact operators, and 
the theory for this kind of operator is much poorer than for the compact 
case. 

Physical restrictions imposed on the non local mass naturally lead us to 
integral equations of the Fredholm type jHj- The 'Fredholm alternative' 
states that, if A = I — A T, where T is a compact operator on a Hilbert space 
H, then the following alternative holds: (a) either Aip = has only the 
trivial solution, in which case Aip = c has a unique solution Vc G H; or (b) 
A<po = has q linearly independent solutions <fi G H. Then A'ip = also 
has q linearly independent solutions (pi G H . In this case Aip = c is solvable 
iff (c,<pi) =0 Vi = l,...,q. 

In the case (b), the general non- homogeneous solution is 

8=1 

where (p p is a particular solution and aj are arbitrary constants. 

On the other hand, when the alternative (a) holds true, this implies that 
the solution shall be unique when c is replaced by <p(0). Solutions correspond- 
ing to c ^ (p(0) are not solutions of the system: (/ — XT)tp — c, c= <p(0), 
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equivalent to the original homogeneous equation (|18|) . and may, therefore, be 
safely discarded. 

Any true solution of the system will also verify a subsidiary equation, 
obtained by setting x 2 = in (JTSJ): 



/+oo 
dy 2 %(0,y 2 )ip(y 2 ) = (24) 
-oo 

for any A ^ 0. This equation shall be true whenever the equations (/ — 
\T)ip — c, c = 99(0) are both true (since it is derived from them), and any 
solution will automatically verify it. 

In our case, we want to study the effect of the non local term, the strength 
of which is controlled by the value of A. Close enough to the local mass 
case, A can be made arbitrarily small, so A -1 is not an eigenvalue of T and 
consequently Aip = has only the trivial solution. Therefore, if T is a 
Fredholm operator, (|18|) will have a unique solution. The functional space 
H to which (p{x) belongs is restricted by the condition 

J dx 2 (^ (x 2 ) ] j (y?(£ 2 )) < 00 , (25) 

because i/j(x) itself has to be normalizable. This becomes a Hilbert space, 
and it contains the zero mode of the local operator (<£> = constant), when 
the scalar product used in H is the one defined by (|23|) . namely, 



(f,g) = J dx 2 (Vo(x 2 )) [f{x 2 )Y g{x 2 ) . (26) 

When Fredholm's hypotheses are satisfied for a general kernel ; jk{x 2 ,y 2 ), 
it is possible to find a perturbative solution. Indeed, writing 

/oo 
dy 2 %(x 2 ,y 2 )(p(y 2 ) , (27) 
-00 

and successively replacing ip(y 2 ) by <^(0) + A J™ 00 dy2%(x2,y 2 )(p(y 2 ) on the 
right hand side, we obtain the Neumann series [H] 

p(0) + A K x {x 2 ) <p(0) H h A n K n (x 2 ) c^(O) + . . . (28) 

where 

/■oo 

K j (x 2 )= / dy 2l W(x 2 ,y 2 ) (29) 
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and 

/oo 
dz 2 %(x 2 ,z 2 )%(z 2 ,y 2 ) 
-oo 

POD 

l (n+1 \x 2 ,y 2 ) = / dz 2 %(x 2 ,z 2 )^ n \z 2 ,y 2 ) . (30) 



It can be shown [§| that the Neumann series converges uniformly to the 
solution <p(x) when 

|A|<-, (31) 

/oo 
dx 2 dy 2 {%(x 2 ,y 2 )f . (32) 
-oo 

So choosing A to satisfy (j31j) . we are allowed to represent (p by the expansion 

<p(x 2 ) = cp{0) + A K x {x 2 ) ip{0) + - ■ ■ + X n K n (x 2 ) ip(0) + ... (33) 

To make the previous analysis more concrete, we consider now an ex- 
ample, based on a specific choice of both the local and non local parts of 
M(x,y) in 0. A natural generalization of the purely local mass case is 
to have m(x 2 ) = Aa(x 2 ) and a jk(x 2 ,y 2 ) which is 'strongly diagonal' and 
symmetric in (x 2 ,y 2 ), i.e.: 

1k{x 2 ,y 2 ) = - a{x 2 ) + a(y 2 ) 5 N (x 2 - y 2 ) , (34) 

where 5n is an approximation of Dirac's delta: 5n(x 2 — 2/2) —* 5{x 2 — y 2 ) 
when iV — > 00. Note that jk(x 2 , y 2 ) — > sign(x2) S(x 2 — y 2 ) when — > 00, so 
that the non local term reduces to a local domain wall mass. 

We adopt 'natural' units such that A = 1, and study the particular case 
\x 2 \ < L = A -1 . Regarding 5n, we use a truncation of the one-dimensional 
completeness relation: 

N 

S N (x 2 ,y 2 ) = Vn(x 2 ) <pi(y 2 ) , (35) 

n=0 

where {f n } is a complete set of functions. We chose ip n {x) to be the harmonic 
oscillator's eigenfunctions. 
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The normalized zero mode corresponding to the local part is 



X2 -, 1 

\X2 



ip {x 2 ) = N exp - I dtm(t) = ^ = e~ |X2 ' , (36) 



vT 



e 



and the corrected zero mode ijjfa) = ^0(^2) ^(^2) is then determined by 
<p(x2). The integral equation for ip{x 2 ) becomes 

ip(x 2 ) = v?(0) + A J dy 2 ; y(x2,y2)ip{y2) (37) 

where 

1 f X2 r 1 

7(^2,2/2) = - / dz 2 e M a{z 2 ) + a(y 2 ) 
* Jo 1 J 

N 

x[X)^)rf(w)]e-'«»l. (38) 

n=0 

From equations ()31j1 and (|3*2*|l . we see that the Neumann series converges, 
in this case, for |A| < 0.0990. So, we assume that A = 0.01, and calculate 
(|37|) perturbatively. Since its expression in terms of analytic functions is not 
very illuminating, we present, in Figure ^ the numerical results of the first 
two iterations, taking (p(0) = 1 and N = 3 . 

In Figure 121 we show the original zero mode (J36|) and the one including 
order-A 2 corrections. We see that the corrected zero mode continues to be 
localized over the domain wall, although it is no longer a symmetric function 
of X2- Besides, this perturbative method introduces only smooth corrections 
in the zero mode, as expected. 

We conclude by summarizing that we have studied a generalization of the 
Callan-Harvey mechanism to the case of a non local mass in 2+1 dimensions, 
showing that for a certain set of assumptions about the non locality, there 
continues to exist a fermionic zero mode. 

We have first considered a quite general non local term, deriving a linear 
integral equation for a function which modulates the zero mode of the local 
case, and accounts for the effect of the non local domain wall. Considering 
a defect of finite size, the Fredholm alternative theorem applies and there is 
a unique, localized chiral zero mode. Moreover, perturbation theory may be 
applied to calculate it: for the concrete example of a strongly diagonal mass, 
we have calculated the first few terms in a perturbative expansion, showing 
that they lead to negligible modifications with respect to the local mass case. 
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Figure 1: f{x2) after two iterations of the Neumann series for the 22-compactificd 
case. The dashed line corresponds to the first iteration. 
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Figure 2: Zero mode profiles for local (dashed line) and non local (full line) masses. 
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